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We report on the dramatic slowing down of the charge carrier dynamics in a quasi-two-dimensional 
organic conductor, which can be reversibly tuned through the Mott metal-insulator transition (MIT). 
At the finite-temperature critical endpoint we observe a divergent increase of the resistance fluc¬ 
tuations accompanied by a drastic shift of spectral weight to low frequencies, demonstrating the 
critical slowing down of the order parameter (doublon density) fluctuations. The slow dynamics 
is accompanied by non-Gaussian fluctuations, indicative of correlated charge carrier dynamics. A 
possible explanation is a glassy freezing of the electronic system as a precursor of the Mott MIT. 


The Mott metal-insulator transition (MIT), where a 
charge gap opens due to electron-electron interactions is 
a key phenomenon in modern condensed-matter physics 
[T], for which the understanding of various fundamen¬ 
tal aspects remains challenging, both theoretically and 
experimentally [ 5 ]. Among them are the nature of the 
anomalous metallic state and (nano-scale) phase sepa¬ 
ration in the vicinity of the Mott transition, the un¬ 
derstanding of the combined effects of electron-electron 
interactions and disorder, and the question of univer¬ 
sality and critical behavior. In recent years, organic 
charge-transfer salts have proven an outstanding class 
of materials with model character for studying these as¬ 
pects of the Mott MIT [JHS]- These materials k-(BEDT- 
TTF)2X, composed of the organic donor molecule 
BEDT-TTF (bis-ethylenedithio-tetrathiafulvalene repre¬ 
senting C6Sg[C2H4]2, abbreviated as ET) and a poly¬ 
meric anion X, are layered systems with a quasi-two- 
dimensional electronic structure and belong to the very 
few examples, where the Mott MIT is purely driven by 
controlling the bandwidth without symmetry breaking 
[ 5 ]. The bandwidth W can be covontrolled either con¬ 
tinuously by applying hydrostatic pressure or in discrete 
steps at ambient pressure by varying the anion X, which 
mimics changes in the ratio W/U, where U is the effec¬ 
tive on-site Coulomb repulsion [ 7 ]. Likewise, replacing 
the eight H-atoms of the ET molecules’ C2H4 ethylene 
endgroups (EEC) in metallic K-(ET)2Cu[N(CN)2]Br {k- 
Hg-Br) by D-atoms (w-Dg-Br) results in a chemically- 
induced shift from the metallic towards the Mott in¬ 
sulating state [ 5 ]. In the T-p or T-X phase diagram 
the MIT is represented by an S-shaped, first-order tran¬ 
sition line which terminates in a second-order critical 
point (pcr./Tcr.), see Fig. 0 Recently, the critical be¬ 
havior of the charge, spin and lattice degrees of freedom 
at (pcr./Tcr.) has been subject of intense research efforts 
[iiMin]. Despite numerous experimental and theoreti¬ 
cal approaches, however, no consensus has been reached 
on the nature of the Mott criticality in «:-(ET)2X and 
the underlying universality class. Furthermore, an in¬ 


vestigation of the critical fluctuations and the dynamic 
properties of the charge carriers at the critical point, in 
particular at low frequencies, is still lacking. From gen¬ 
eral arguments, a critical slowing down of the fluctuations 
is predicted for a second-order phase transition [IB] but 
has not yet been observed at the Mott critical endpoint 
0 - 

In this letter, we present a systematic study of the dy¬ 
namical properties of the electrons at the Mott transi¬ 
tion and critical endpoint using fluctuation spectroscopy, 
a technique, where the dynamics of charge carriers is 
studied without injecting additional electrons into the 
system. For our study, we have chosen the material n- 
[(Hg-ET)o.2(Dg-ET)o.g]2Cu[N(CN)2]Br (K-Dg/Hg-Br in 
short), where the partial substitution of the ET molecules 
with their deuterated analogues places the sample on the 
high-pressure (metallic) side of the phase diagram very 
close to the critical pressure Pcr. of the MIT m Upon 
fine-tuning the material through the transition we find an 
electronic-correlation-induced enhancement of the low- 
frequency resistance fluctuations, which strongly increase 
and nearly diverge when approaching the critical point, 
accompanied by a dramatic slowing down of the fluctua¬ 
tion dynamics. At the temperatures, where the fluctua¬ 
tions are largest, non-Gaussian fluctuations are observed. 

«:-[(H8-ET)o.2(Dg-ET)o.g]2Cu[N(CN)2]Br has been 
grown as single crystals by electrochemical crystal¬ 
lization [TH]. Low-frequency fluctuation spectroscopy 
measurements have been performed in a five-terminal 
setup using a bridge-circuit ac technique m described 
in detail in [20] . Resistance R and resistance noise power 
spectral density (PSD) S'i?(/) = 2(|JU(/)p)//^ |2T| have 
been measured perpendicular to the conducting layers, 
where I is the applied current. For all measurements 
discussed here, we have observed excess noise of general 
Sr oc I//“-type (see |22l (23) for typical spectra), 
characterizing the intrinsic resistance (conductance) 
fluctuations of the sample. 

Figure[^a) shows the resistance of K-Dg/Hg-Br for var¬ 
ious rates of cooling q = dT /dt through the temperature 
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FIG. 1. (Color online) (a) Resistance measurements of k- 
Dg/Hg-Br for various cooling rates q. (b) Arrow indicates the 
position of the slowly-cooled, pristine sample in the generic 
phase diagram. Open circles indicate values taken from the 
literature of the second-order critical end point of the first- 
order Mott MIT (thick solid curve), see |I7) and references 
therein. The change of the sample position with q is schemat¬ 
ically indicated by the colored lines. Blue represents fully 
deuterated K-Dg-Br |23| . 


Tg « 75 — 80K of the glass-like EEG ordering [211 E5] . 
which results in a frozen non-equilibrium occupation of 
the two possible EEG conformations [21]. The trans¬ 
formation from metallic to insulating behavior with in¬ 
creasing q{Tg), see Fig.[^ has been widely used in K-Br 
compounds with different type and degree of deuteration 
of the ET molecules [I71IIS1I271I2H]- The effect originates 
in an anisotropic change of the in-plane lattice parame¬ 
ters at Tg, which changes the relevant transfer integrals 
such that more rapid cooling leads to slightly smaller 
bandwidth (and W/U ratio) [T71 [21]. In Fig.j^a), the 
sample which has been annealed for 5 h in the vicinity 
of Tg, shows essentially metallic behavior. In contrast, 
the resistance curves taken during continuous cool-down 
measurements with increasing q{Tg) reveal the typical 
re-entrant behavior, where the S-shaped first-order Mott 
transition line is crossed twice. In between the slowest 
and fastest cooling-rate of the experiment, the sample has 
crossed the second-order critical endpoint {pcr ., Ter .) of 
the Mott transition. All cool-down measurements shown 
here reveal a superconducting transition at low temper¬ 
atures, which is of inhomogeneous, percolating type, due 
to the coexistence of metallic/superconducting and in¬ 
sulating/antiferromagnetic phases [1111211 ISO] • We note 
that this fine-tuning of the electronic bandwidth by con¬ 
trolling the cooling rate is reversible, i.e. upon warming 
above Tg the frozen EEG glass melts and the lattice re¬ 
laxes. 

After the continuous cool-down runs with different 
rates q, noise measurements have been performed in 
each case upon warming the sample in discrete tem¬ 
perature steps. Each data point in Fig.[^ represents 
one 1 //“-type spectrum in the frequency range between 
1 mHz and 100 Hz, where the magnitude of the fluctu¬ 


ations, Sfi/R^^f = I Hz), and the frequency exponent, 
a = dlnSa/dlnf, have been evaluated. Figurej^a) - (c) 
shows the resistance noise PSD for K-Dg/Hg-Br in com¬ 
parison to the fully deuterated compound K-Dg-Br (data 
from [21] ), which is located slightly on the insulating side 
of Mott transition, see Fig. Bb). The data are normalized 
to the broad maximum in Sr/R^ at around 100 K, which 
is due to the coupling of the charge carriers to the fluc¬ 
tuating EEG of the ET molecules [211HH [H] , which un¬ 
dergo a glass-like ordering transition at Tg ^ 75 — 80 K as 
described above. Since this feature is of structural origin, 
it is independent of the samples’ position in the gener¬ 
alized phase diagram. Figure]^ a) compares the slowly- 
cooled fully and partially deuterated compounds K-Dg- 
Br and K-Dg/Hg-Br, respectively, showing essentially the 
same behavior at high temperatures T ^ Tg. For k- 
Dg-Br, the arrow indicates a sharp peak in the noise at 
T = 36K (however, small in comparison to the broad 
maximum at 100 K), which grows larger at frequencies 
smaller than IHz (see Figs. 1 and 3 in [23] )■ This fea¬ 
ture, which is absent in the fully-hydrogenated metal¬ 
lic system K-Hg-Br (not shown), has been interpreted as 
a sudden slowing down of the charge carrier dynamics 
due to the vicinity of the Mott critical endpoint [23]. In 
contrast to K-Dg-Br, we observe a larger overall noise 
level for K-Dg/Hg-Br and even a slight increase below 
Tg 75 — 80K, i.e. where the conformational vibrations 
of the EEG are frozen out. This can be understood con¬ 
sidering the sample’s position in the phase diagram closer 
to the Mott insulating phase resulting in a correlation- 
induced increase of the low-frequency fluctuations, since 
the charge carriers already tend to localize, an effect, 
which possibly is enhanced by weak static disorder m- 
Likewise, electronic phase separation in this region may 
enhance the low-frequency fluctuations [30l [32] . Most 
striking, however, is the sharp and - compared to K-Dg- 
Br - strongly enhanced noise peak at T ^ 33K for the 
K-Dg/Hg-Br sample cooled with q = 0.5K/min. 

These large low-frequency fluctuations become en¬ 
hanced even more dramatically for increasing cooling 
rates oi q = 2K/min and 5K/min shown in Figs.j^b) 
and (c), respectively [33]. Note the rescaled ordinate in 
Figs.j^c) demonstrating the very large noise level ob¬ 
served for the data taken at 5K/min, i.e. this particular 
position in the phase diagram, while for q = lOK/min 
the noise peak decreases again. This observation (i) rules 
out a mere effect of increasing lattice disorder [26], and 
(ii) suggests that the curve q = 5K/min is the one clos¬ 
est to the critical endpoint of the Mott transition, with 
the highest noise level at T^. = 33 K. The frequency ex¬ 
ponent a{T) for this curve, shown in Fig.j^d), increases 
upon approaching Ter. from values of 0.8 at T = 48K to 
a = 1.6 at 31K, corresponding to a strong shift of spec¬ 
tral weight to low frequencies. In general, large values 
of a > 1.4 are typical for a system far from equilibrium, 
as we will discuss below. At the lowest frequency of our 
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FIG. 2. (Color online) (a) - (c) Resistance noise (PSD) taken at / = 1 Hz of rt-Ds/Hg-Br compared to K-Dg-Br (slowly cooled, 
same data as in m ) normalized to the value of the local maximum around 100 K for increasing cooling rates (a) q — 0.5, (b) 2, 
(c) 5 and lOK/min. Note the rescaling of the ordinate in (c). (d) Frequency exponent a vs. T of the l//“-type noise spectra 
for q = 5K/min. Line is a description with a model of non-exponetial kinetics (DDH model) as described in the text, (e) 
Relative noise level an = f x Sr/ B? vs. f vs. T for the same cooling rate. 


experiment of / = 1 mHz, the noise level increases by 
more than three orders of magnitude upon approaching 
the critical point from above. In the study of dynamical 
critical phenomena the existence of such large times is 
known as critical slowing down of the order parameter 
relaxation rate [H], where the length scale which mea¬ 
sures the correlations of the order parameter becomes 
infinite upon approaching a second-order phase transi¬ 
tion. The time evolution of the order parameter becomes 
very slow since the correlated regions in the sample in¬ 
volve very many degrees of freedom. Therefore, the ki¬ 
netics of the order parameter near a critical point causes 
a slow dynamics of a macroscopic system, in our case 
manifested in the diverging resistance/conductance fluc¬ 
tuations. The expected singular nature of the critical 
point as probed by order-parameter fluctuations is high¬ 
lighted in Fig.j^e), where we plot the relative noise level 
a]i{f,T) = / X Sr/B?, a dimensionless number char¬ 
acterizing the strength of the fluctuations. Our results 
demonstrate the expected divergence of the doublon (two 
electrons sitting on the same atomic orbitals) suscepti¬ 
bility and the critical slowing down of doublon density 
fluctuations 

Besides the strongly increasing amplitude of the low- 
frequency fluctuations at the critical point, another im¬ 
portant aspect related to the dynamic criticality can be 
deduced from our data. As we have shown in our ear¬ 
lier studies [551 US] , the 1 //-type fluctuation properties 
of the organic charge-transfer salts k-(ET )2 can be very 


well described by a model of non-exponential kinetics 
[36] ■ where the l//“-noise originates in the superposi¬ 
tion of a large number of independent, i.e. Gaussian, so- 
called ’fluctuators’. Within this model, first introduced 
by Dutta, Dimon and Horn (DDH) [37], a distribution of 
activation energies D{E) of the individual fluctuations 
determines the observed frequency and temperature de¬ 
pendence SR{f,T). The assumptions of the model can 
be checked for consistency by comparing 


aDDH(E) = 1 — 


1 

In 27r/To 


( d\nS{f,T) 

51nT 



( 1 ) 


with the measured data a{T) [5B]. The blue line in 
Fig.j^d) shows aDDH('E), which describes the measured 
data very well, except for an obvious deviation just above 
the critical point at T^r., where the noise level peaks. We 
note that such a deviation, which we also see for the 
q = lOK/min curve, is not observed for K-Dg-Br |23j 
and for the smaller cooling rates of the present sample 
K-Dg/Hg-Br (not shown). An obvious interpretation is 
that the assumption of independent two-level fluctuators 
is simply not valid in a narrow temperature range around 
Ter., which is a first hint to non-Gaussian fluctuations 
near the critical point. 

In order to further explore this matter, we introduce 
the so-called second spectrum f, T), which is the 

power spectrum of the fluctuations of Sr//, T) with time, 
i.e., the Fourier transform of the autocorrelation func¬ 
tion of the time series of Sr^J) [35], introducing an addi- 
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tional frequency /2 related to the time over which Sr^/) 
fluctuates [39]. probes a fourth-order noise 

statistics and therefore deviations from Gaussian behav¬ 
ior: it is independent of the frequency /2 (’’white”) if the 
fluctuations are uncorrelated, i.e. caused by independent 
two-level systems. In contrast, oc 1//^ for correlated 
(interacting) fluctuators (301 SI]- As shown in Fig.for 
K-Dg/Hg-Br cooled-down with q = lOK/min, i.e., even 
slightly beyond the critical point, the second spectrum 
is ’white’ above and below the temperature of the noise 
peak in the first spectrum. At T = 34K, where the noise 
peaks, however, is frequency-dependent with a large 
exponent = 1.35, see black line in Fig.|^ Clearly, close 
to the critical point, where the fluctuation dynamics is 
slow, this spectral wandering indicates correlated fluctu¬ 
ations, which rules out independent two-level systems as 
the source of the noise. 
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FIG. 3. (Color online) Second spectra S^^\f 2 ,f) vs. f 2 /f 
of K-Dg/Hg-Br for q = lOK/min at three different tempera¬ 
tures. Horizontal lines are guides for the eyes. At T = 34 K, 
where the first noise spectra SR/R^{f) peaks at / = 1 Hz, 
see Fig.j^c), the data at low frequencies can be fitted to 
oc I //2 yielding /3 = 1.35 (black line). At all temper¬ 
atures the spectra measured at different frequency octaves 
are scale invariant, i.e. do not depend on the scale of /. 

In general, non-Gaussian fluctuations caused by slow 
modes corresponding to large scales are expected near 
a critical point of a second-order phase transition and 
have been observed, e.g., in liquid crystals jUj- However, 
since non-ergodic dynamics is found in a wider range of 
bandwidths around the critical region, i.e., for q = 5, 
lOK/min and faster cooling rates, it is interesting to 
compare our findings to noise measurements of other sys¬ 
tems showing metal-insulator transitions. For a 2D MIT 
in Si inversion layers, a similar dramatic slowing down 
of the electron dynamics is observed |43j . which is in¬ 


terpreted as a glassy freezing of the electron system, in 
agreement with the onset of non-Gaussian fluctuations 
already in the metallic phase. Similar observations for 
the 3D Anderson-Mott MIT in P-doped Si are likewise 
interpreted as a glassy freezing of the electrons preced¬ 
ing their localization |44j . and the intriguing possibility 
that such correlated dynamics may be a universal feature 
of MITs, irrespective of the systems’ dimensionality, is 
raised. We note that in these systems the noise diverges 
only for T —> 0, whereas for the present system a drastic 
increase and slowing down of the charge-carrying fluctua¬ 
tions is observed for the first time at a finite-temperature 
critical endpoint. A glassy freezing of the charge carriers 
near the critical point may be a possible scenario also for 
the present strongly-correlated materials US] [IS]- Here, 
a residual structural disorder potential and/or the inher¬ 
ent frustration of the triangular lattice geometry of the 
FT molecules forming dimers could result in a large num¬ 
ber of metastable states, which the correlated electrons 
have to overcome. General models for such glassy dynam¬ 
ics applied to spin glasses are often invoked also for the 
charge degrees of freedom m- In a model of interacting 
droplets/clusters f) should be a decreasing func¬ 

tion of / for constant / 2 // |40l[47]. The observed scale 
invariance shown in Fig.[^ (note the different colors for 
each frequency octave), however, points to a hierarchical 
picture, where the system wanders collectively between 
metastable states related by a kinetic hierarchy [ID]. The 
large exponent of /3 = 1.35 indicates that only few such 
states are visited during the time of measurement. From 
the theoretical point of view, a Coulomb glass behavior 
in the vicinity of the Mott transition remains controver¬ 
sial. For example, a self-generated glass transition caused 
by the frustrated nature of the interactions (and not re¬ 
lated to the presence of quenched disorder) is predicted in 
doped Mott insulators |3S]. In [35], however, the metallic 
glassy phase is suggested to be suppressed for Mott local¬ 
ization and to become stabilized only for increasing dis¬ 
order, in apparent agreement with experimental results 
on Si inversion layers with different degrees of disorder. 

Our method applied to the present highly-tunable 
molecular conductors opens the door to systematically 
explore the dynamic critical properties of strongly- 
correlated electrons in the presence of disorder and 
frustration. Whether a universal electron glass state 
emerges in the vicinity of the finite-temperature critical 
endpoint of the Mott transition in k-(ET) 2 X or non- 
ergodic dynamics and the relevance of higher statistical 
moments is an inherent characteristics of the proximity 
to the critical point of the second-order phase transition 
will be seen in future systematic studies of fluctuation 
spectroscopy for different bandwidths and disorder 
levels. Furthermore, in view of the large anomalies in 
the sound velocity at the critical point, which is consis¬ 
tent with a diverging compressibility of the electronic 
degrees of freedom [3H], and in the lattice response at 
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the critical point [JS], the coupling of the diverging 
low-frequency fluctuations to the lattice degrees of free¬ 
dom are an important field of investigations in the future. 
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